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Abstract 

Vortons are closed loops of superconducting strings carrying current and charge. A 
formalism has been developed to study vortons in terms of an elastic string approxi- 
mation, but its implementation requires knowledge of the unknown equation of state, 
relating the string tension to the energy per unit length. Recently, a planar analogue 
of the vorton, known as a kinky vorton, has been introduced. In this paper we derive 
an exact formula for the equation of state of a kinky vorton and use it to calculate 
the properties of the associated elastic string, such as the transverse and longitudinal 
propagation speeds. In particular, the elastic string approximation predicts a compli- 
cated and highly non-trivial pattern of intervals of instability, which we are able to 
confirm using full field simulations. The implications of the results for vortons are also 
discussed. 



1 Introduction 



Superconducting cosmic strings were introduced by Witten [13], who realized that if the 
field of a cosmic string is coupled to another complex scalar field then a non-dissipative cur- 
rent can fiow along the string. A vorton [8j is a closed loop of superconducting cosmic string 
that carries both current and charge, which may provide a force to balance the string tension 
and prevent its collapse. As vortons have a number of possible cosmological consequences 
it is of considerable interest to determine their stationary and dynamical properties. 
Recently, results have been presented [2] on the first numerical construction of vortons 
in the simplest theory, namely, the global version of Witten's U{1) x U{1) theory. These 
results demonstrate that a range of stationary circular vortons exist, that are stable to axially 
symmetric perturbations. However, unstable modes associated with non-axial perturbations 
were found to exist. 

A formalism has been developed [H [6] to study vortons and their stability in terms of 
an elastic string approximation, but its implementation requires knowledge of the unknown 
equation of state, relating the string tension to the energy per unit length. Attempts have 
been made pjj to compute the equation of state numerically from field theory simulations 
using particular values of the field theory parameters, but obviously this is far from ideal. 

The equation of state can be derived from the action density on the string worldsheet, 
but one needs to know how this depends on the quantity x = tu^ — A;^, where uj and k are the 
frequency and twist rate of the phase of the condensate field. Witten originally suggested [13] 
that the usual Nambu Goto action density, which is simply a constant, could be replaced by a 
form linear in x- However, this linear form produces results which are qualitatively incorrect, 
even for small X; since the propagation speed of longitudinal perturbations depends on the 
second derivative of the action density with respect to x- This prompted suggestions [3 O [9] 
for alternative nonlinear forms for the assumed equation of state. 

Recently, a planar analogue of the vorton, known as a kinky vorton, has been introduced 
and shown to possess many of the features expected of a vorton [T]. A significant advantage 
of the kinky vorton is that several approximations required in the study of vortons can be 
replaced by exact results. In this paper we exploit this advantage to derive an exact formula 
for the kinky vorton equation of state, and make a comparison with the previously proposed 
approximations to the vorton equation of state. The exact equation of state is then used 
to calculate properties, such as the transverse and longitudinal propagation speeds, of the 
elastic string description of the kinky vorton. In particular, the elastic string model predicts 
a complicated and highly non-trivial pattern of intervals of instability, which we are able 
to confirm using full field simulations. These analytic results provide an understanding of 
similar instabilities found for vortons using numerical simulations of the full field theory 
dynamics [2]. 
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2 The action density on the string worldsheet 



The kinky vorton Lagrangian density in (2+l)-dimensions is given by [Ij 

c = + d.ad^a - ^(0^ - vir - ^(i-p - vir - + ^vt (2.1) 

where and a are real and complex scalar fields respectively, with i]^,-)]^, X^p, X„, f3 all real 
positive constants. 

This theory can be obtained from the global version of Witten's U{1) x U{1) theory [13j 
by a trivial dimensional reduction from (3+l)-dimensions to (2+l)-dimensions, followed by 
a restriction that one of the complex scalar fields is real. 

The theory has a global Z2 x f/(l) symmetry and the parameters of the model can be 
arranged so that in the vacuum the Z2 symmetry is broken, = ±r/<^ 7^ 0, while the U{1) 
symmetry remains unbroken, \a\ = 0. For this symmetry breaking pattern there exist kink 
strings constructed from the field. If the infinite kink string lies along the y-axis, then it 
is given by the solution 



= r/,tanh(^!^^^), a = 0. (2.2) 

The situation of interest is when a condensate of the a field carrying current and charge 
forms in the core of the kink string. For the infinite string given above, such a condensate 
field takes the form 

a = e*('^*+^^Vl, (2-3) 
where \a\ is a function of x only with |cr|— s>Oas|x|-^oo. The constant k describes the rate 
of twisting of the condensate along the string, though in the literature this is referred to as 
winding, rather than twisting, so we will stick to this common convention. 

A non-zero value of u induces a charge Q associated with the global U{1) symmetry, and 
the winding k generates a current along the string. It is easy to see that charge and current 
will have opposite effects on the string, so it is useful to introduce the combination 

X = oo^- k\ (2.4) 

In the literature solutions with x = ^I's termed chiral, whereas solutions with x > ^-I's 
referred to as electric and those with x < are called magnetic [H [10] . 
To obtain exact solutions for the whole range of x "we set [1] 

2j3 = \^ = K = A. (2.5) 

With this choice there are three remaining parameters, which may be chosen to be t]^ and 
the two combinations 

2 _ 2 _ [Va 



In fact, the parameters and merely set energy and length units and can therefore 
be set to unity without loss of generality. To restore these parameters in the following the 
scalings L 1— > m^rj^L and x ^ xl'^,^ need to be applied. 
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The only significant parameter that remains is a, which must he in the interval a G (|, 1) 
if chiral solutions are to exist. The range of x is 

X- = ^- "<X<l-« = X+, (2.7) 
and the associated exact solutions are 



= tanh(2;v/l - a-x), a = e^('^*+'^^V2a - 1 + 2xsech(xv/l - « - x)- (2.8) 

Substituting the solution (12.81) into the Lagrangian density (12.11) and integrating over x 
yields the result 

'(2a + l)(a-l) + x(2x + 4a-l)| (2.9) 



3^1 - a - X 

for the action density on the string worldsheet. 

This should be contrasted vith various proposed forms, such as the linear approximation 

Ly,^ = -vr? ^ Kx. (2.10) 

due to Witten [13] , and the logarithmic approximation [7J [9] 

= -vr? - \ml\og{\ - 8lx). (2.11) 

The parameters vn? and k in the linear approximation (12.101) can be determined from the 
exact expression (12.91) by matching the function and the derivative at x = 0. This gives 

9 4 , , , 4a - 2 , , 

= -{2a + l)Vl - a, = , (2.12) 

3 -y/l — a 

One approach to fixing the additional parameters ml and 6^ in the logarithmic approximation 
is by matching the derivative at x = 0, and requiring the singularity in (12. lip to occur at 
the extreme electric limit x = X+) which gives 



m 



2 



4(2a-l)v^r^, 6l = l/{l-a). (2.13) 



In Figure [T] the exact result (12. 9p is compared with the linear (I2.10p and logarithmic 
(12. lip approximations for the parameter value a = 3/4, with other parameter values pro- 
ducing similar results. This shows that the linear approximation is poor, but the logarithmic 
approximation is reasonably good, and therefore should produce results which are, at least, 
qualitatively correct. 



4 







_l 



-1.5 - 



-0.5 - 



-1 - 




-2 



-0.25 -0.2 -0.15 -0.1 -0.05 







0.05 



0.1 



0.15 



0.2 



0.25 



X 



Figure 1: The exact action density on the worldsheet L (sohd curve) as a function of Xi to- 
gether with the hnear approximation Liin (dashed curve) and the logarithmic approximation 
Liog (dotted curve). This is in the theory with a = 3/4. 

3 The equation of state 

From the Lagrangian (12.91) the string tension T and energy per unit length U can be 
calculated as follows J4j. 

The energy-momentum tensor associated with the Lagrangian density (12.11) is 



Integration over the string cross-section gives the macroscopic tensor 



where a, 6 G {t, y}. The energy per unit length U and the string tension T are the eigenvalues 
of T"'', and are given hj U = T** and T = — T^^, in a frame in which T"^ is diagonal. 
Applying the above procedure to the Lagrangian density (12. ip produces 




(3.2) 



T" = 2uj^E2 - L, 



(3.3) 



where 




(3.4) 
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The frame in which T"'' is diagonal is obtained by setting A; = 0, if x > 0, or by setting 
= if X < 0. 

In the electric regime, that is x > 0, setting k = implies that x = and then (13.31) 
becomes 

U = 2x^2 ~L, T = -L. (3.5) 

Conversely, in the magnetic regime, that is x < 0, then setting u = means that x = 
and hence (13. 3p gives 

U = -L, T = 2x^2 - L. (3.6) 

Of course, both forms agree in the chiral limit x = where U = T = —L. 

Using the explicit expressions (12.91) and (13.41) for L and S2 the above equations become 

(2a + !)(!-«)- x(2x + 4a - 1)) /^/1-a-x if X > 

T={ } { (3.7) 

^ ^ (2a + 1)(1 - a) + 2x(2x + a - 1) ] /^1-a-x if X < 

and 

' (2a + 1)(1 - a) + 2x(2x + a - 1) ) / ^1 - « - X if X > 

U={ } { (3.^ 

(2a + 1)(1 - a) - x(2x + 4a - 1) j /^l - « - X if X < 0. 

Using these formulae it is easy to calculate the behaviour of the tension and energy in the 
three limits x = (chiral), x ^ X+ (extreme electric) and x ~^ X- (extreme magnetic). 

If X = then T = U = = |(2a + l)^/T^. As x ^ X+ then T ^ and f/ ^ 00. 
Finally, as x X- then T — ^ 4-\/2/3 and U 4-\/2/3, so again the tension and energy are 
equal, and independent of a. 

Figure [2] displays the tension T (solid curve) and the energy per unit length U (dashed 
curve) as a function of x, for the theory with a = 3/4. This value of a has been employed in 
previous work since the range of x is then symmetric around zero, explicitly x+ = 1/4 and 
X- = —1/4, allowing a range of electric and magnetic kinky vortons. Plots for other values 
of a share the same qualitative features. All graphs in this paper are plotted for the theory 
with a = 3/4, but again the qualitative features are independent of a. 

Note that the tension is never negative, which agrees with numerical computations in 
(3+l)-dimensions, confirming that there are no spring states [I2], despite earlier claims in 
the literature (see the discussion and references in [12j). 

The equation of state is the relation between T and U, and this is given implicitly by 
combining the formulae (13.71) and (13.81) . Figure [3] presents the equation of state in graphical 
form, that function of U. For comparison the upper dashed line is the Nambu- 

Goto equation of state T = U, and the lower dotted curve is the self-dual equation of state 
[5] T = m'^/U, with the characteristic property that the equation of state is identical in the 
magnetic and electric regimes. The self-dual model appears to be a reasonable approximation 
in the electric regime, but as discussed in [7] and demonstrated below, the self-dual model 
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Figure 2: The tension T (solid curve) and the energy per unit length U (dashed curve) as a 
function of for the theory with a = 3/4. 



does not capture the important qualitative features associated with propagation speeds, even 
in the electric regime. 

Figure [3] has a remarkable similarity to the graphical equation of state computed nu- 
merically for superconducting cosmic strings in (3-|-l)-dimensions [TT]. This suggests that 
our explicit exact formulae derived in (2+l)-dimensions also provide a good description of 
the (3-|-l)-dimensional system, where results are only available numerically. Note that the 
most interesting dynamics of a string loop takes place in the plane of the loop, therefore 
the reduction to (2+l)-dimensions is still expected to capture the most important degrees 
of freedom. 

String dynamics depends crucially on the propagation speeds of transverse and longitudi- 
nal (sound-like) perturbations. The equation of state allows the calculation of the transverse 
speed ct and the longitudinal speed cl via the formulae 

Ct = jj^ = (3-9) 

Substitution of the tension and energy expressions (13. 7p and (13. 8p into the speed formulae 
dSSD gives 



(3.10) 



2sign(x) _ (2a + 1)(1 - a) - x(2x + 4a - 1) 
(2« + l)(l-a) + 2x(2x + «-l)' 

2sign{x) ^ (2a + 2x - 1)(1 - g - x) /o 
^ (2«-l)(l-«)-2x(2x + 3a-3)' ^ ' ^ 
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Figure 3: The tension T as a function of the energy per unit length U (sohd curve), for 
a = 3/4. The upper portion of the curve is the magnetic regime (x < 0) and the lower 
portion is the electric regime (x > 0). For comparison the upper dashed line is the Nambu- 
Goto equation of state T = U, and the lower dotted curve is the self-dual equation of state 
T = m^/U. 



with cr = cl = 1 if X = 0. 

The string is unstable if either of the speeds are imaginary, that is, if either cf^ < or 
c\ < 0. It is easy to see that cf, > for all x ^ (X-? X+)^ ^'iid ct — >■ as x — X+- Similarly, 
it is easy to show that c| > for all x ^ [0,X+)? a^id again — > as x ^ X+- However, 
c| < if X < Xci where the critical value for the onset of instability is given by the vanishing 
of the denominator in fl3.1ip . leading to 

Xc = ^(1 -a)- ^V(l-a)(5-«). (3.12) 

For the parameter value a = 3/4 this gives Xc = (3 — vT7)/16 ~ —0.07. This is in excellent 
agreement with numerical field theory simulations, performed for the value a = 3/4, which 
estimated the critical value to be —0.08 [1]. 

In Figure m the speeds ct (solid curve) and cl (dashed curve) are plotted as a function 
of X; for C( = 3/4. Plots with other values of a are qualitatively similar. 

In the chiral limit x = 0? ^''^^ extreme electric limit x X+) the string is transonic, 
that is, Ct = c^. For all other values of x the string is supersonic, that is, ct > cl- Such super- 
sonic behaviour has been observed in numerical computations based on (3+l)-dimensional 
field theories [TT]. 
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Figure 4: The transverse speed ct (solid curve) and the longitudinal speed cl (dashed curve) 
as a function of Xi for the theory with a = 3/4. 



The simple hnear approximation to the action (12.101) is quahtatively incorrect, since it 
predicts that the string is everywhere subsonic, that is, ct < cl, which we have seen is never 
the case. Models based on a self-dual equation of state predict that the string is everywhere 
transonic, which is also false. These failures motivated attempts to construct improved 
actions [3 El [9] , of which the logarithmic approximation (12.111) is an example. As we have 
seen, the exact action for kinky vortons has a simple form (I2.9p . but an action of this type 
does not appear to have been suggested previously. Given the remarkable similarity between 
the results obtained using (I2.9P and numerical computations in (3+l)-dimensions, it seems 
likely that an action of the form (I2.9P will provide a good description for vortons. 

4 Intervals of instability 

For a circular elastic string of radius R, the frequencies Qn, of linear perturbations with 
fourier mode n, can been related to the transverse and longitudinal speeds [B]. Explicitly, 
the scaled frequency Un = ^uR/ct satisfies the cubic equation 



a3^n + <^2l^l + O-l^n + ^0 = 0, 



(4.1) 
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where 



ao = 2(4 -4)(n2 - l)n (4.2) 
ai = 44(1 - 4) (n^ - 1) - (1 + 4) (4 - 4) (n^ + 1) (4.3) 

aa = 24(4 - 4 - 2(1 - 44))^ (4-4) 
ag = 4(l + 4)(l-44). (4.5) 

The derivation of this formula in [6] assumes that the string moves in a three-dimensional 
space, but the result applies equally well to a string in two-dimensional space because the 
modes perpendicular to the plane of the circular string decouple and are irrelevant for the 
stability analysis. 

Instability is characterized by a complex root of the cubic (14.11) . It is easy to show [6] that 
all roots are real for n = and n = 1, simply as a consequence of the causality restrictions 
< 4 — 1 < 4 ^ 1 non-zero speeds. 

For axially symmetric perturbations (n = 0) the solutions of the cubic are a trvial zero 
mode and 

2 ^ 24(i-4) + (4 + 4)(i-4) 

° 4(i + 4)(i-44) ■ ^ ' 

The explicit speed formulae ( 13.10p and ( 13. lip allow this frequency to be calculated for any 
kinky vorton given the values of x and R. As an example, a kinky vorton is reported in [1] 
with Q = 1500 and N = 84, producing the quantities x = 0.0472 and R = 154.6. Using 
these values in the above formulae gives a frequency Qq = 0.00644. This can be compared 
with the results of full field dynamics, using the numerical approach described in detail in 
[1]. Perturbing this kinky vorton by an axially symmetric mode that increases the radius 
by 1% produces an oscillation of the radius. A spectral analysis of this oscillation yields the 
frequency 0.00658, which is very close to the above predicted value for Qq. 

For a given mode n, the critical values of x niarking the limits between stability and 
instability correspond to the values at which the cubic (14.11) has a repeated root. This is 
given by the vanishing of the resultant 

Res(a3Z/^+a2i^„+ai^'„-|-ao, 3a3Z/„+2a2^'n+ai) = 27ar^aQ — 18a^aoaia2-~asa2a^+Aa2a3aQ+Aa^a^. 

(4.7) 

In Figure [5] we indicate the x intervals of instability for modes up to ri = 10, for the 
theory with a = 3/4. Note that the range of x for the plot runs from the value at which 
the longitudinal speed vanishes Xc = (3 — vT7)/16 (in which case there is automatically 
instability). 

This plot reveals a complicated and highly non-trivial structure for the instability inter- 
vals as a function of mode number. In particular, if a kinky vorton is to be stable to all 
perturbation modes then there is an extremely limited range in the magnetic regime (x < 0). 

Note that the full field dynamical simulations of kinky vortons presented in [1] included 
the perturbation of a solution with x = —0.04 that did not decay. However, in this case 
the perturbation was generated by the square boundary and therefore corresponds to modes 
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Figure 5: The x intervals of instability for perturbations with fourier modes up to n = 10, 
for the theory with a = 3/4. In this plot the lower limit of x is equal to the value at which 
the longitudinal speed vanishes. 



that are a multiple of 4. This result is consistent with Figure [5] since the unstable modes for 
X = —0.04 are n = 5 and n = 6. 

Most electric vortons (x > 0) are stable to the instabilities discussed above and the 
results are consistent with full field simulations pj. The elastic string model predicts an 
instability for a small range of electric vortons around x = 0.05, under a perturbation with 
mode n = 2. A suitable kinky vorton solution to investigate this issue is the one discussed 
above with Q = 1500 and = 84 since this has x = 0.0472 and R = 154.6. The elastic 
string model predicts an instability to n = 2 perturbations with the associated frequency 
Q2 = 0.0064 + 20.00023. The period of kinky vorton oscillations is already large, but the 
imaginary part of Q2 is a further order of magnitude smaller than the real part. This implies 
that the timescale over which the instability manifests itself is very large and long simulations 
are required to observe the instability. 

The numerical approach described in [T] has been applied to study the stability of the 
above kinky vorton. The perturbation consists of a squashing by 2% along the x-axis and a 
stretching by the same factor along the y-axis, in order to preserve the total charge Q. This 
elliptic deformation induces an n = 2 perturbation and the resulting evolution is presented 
in Figure [6], where the solid curve represents the radius of the kinky vorton, as measured 
along the x-axis. It is clear that the perturbation grows with time, confirming that the kinky 
vorton is unstable to modes with n = 2. The dashed curve displays the predicted growth rate 
of the envelope, that is, it is a curve with the predicted growth exp(tS5(fi2))- Once again. 
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Figure 6: The evolution of the radius, Rx (sohd curve) calculated along the x-axis, under 
an n = 2 perturbation, for a kinky vorton in a regime where instability is predicted. For 
comparison, the expected growth of the envelope is also presented (dashed curve). 



this shows an excellent agreement with the elastic string analysis. 

The results described above, together with other similar simulations, including those 
presented previously in [1], confirm the validity of the elastic string description. The gen- 
eral conclusions are that most magnetic kinky vortons are unstable to generic non-axial 
perturbations, whereas the stability of electric kinky vortons has a crucial and complicated 
dependence upon the parameters of the kinky vorton. 

5 Conclusion 

In this paper we have derived an exact formula for the action density on the string word- 
sheet in the kinky vorton model. Using this result we have determined the equation of state 
and hence explicit expressions for the transverse and longitudinal propagation speeds in the 
elastic string description of a kinky vorton. This has allowed an explicit analytic study of 
the stability of kinky vortons, revealing a highly non-trivial pattern of intervals of instabil- 
ity. The analytic results are in excellent agreement with numerical results from full field 
simulations, confirming the validity of the elastic string description. 

The exact kinky vorton results share the same qualitative features found in numerical 
computations of vortons, which suggests that the kinky vorton form of the action density on 
the string worldsheet should provide a good approximate description for vortons. 



12 



A recent numerical investigation of vortons [2] found instabilities to non-axial pertur- 
bations which are very similar to the instabilities described here for kinky vortons. The 
kinky vorton results show that the existence of an instability is not a generic feature, but 
rather has a crucial and non-trivial dependence on the properties of a particular kinky vor- 
ton. Numerical simulations of vortons requires considerable computational resources and 
unfortunately this restricts investigations to quite a limited region of parameter space. The 
results of the present paper suggest that the instabilities found in vorton simulations [2] 
may not exist for all vortons. Certainly, an important observation is that results found in a 
limited region of parameter space are unlikely to be generic throughout the parameter space 
of vorton solutions. 
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